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ABSTRACT 

The formation of jets such as dynamic fibrils, mottles, and spicules in the solar chromosphere is one 
of the most important, but also most poorly understood, phenomena of the Sun's magnetized outer 
atmosphere. We use extremely high-resolution observations from the Swedish 1-m Solar Telescope 
combined with advanced numerical modeling to show that in active regions these jets are a natural 
consequence of upwardly propagating slow mode magnetoacoustic shocks. These shocks form when 
waves generated by convective flows and global p-mode oscillations in the lower lying photosphere 
leak upward into the magnetized chromosphere. We find excellent agreement between observed and 
simulated jet velocities, decelerations, lifetimes and lengths. Our findings suggest that previous obser- 
vations of quiet sun spicules and mottles may also be interpreted in light of a shock driven mechanism. 
Subject headings: magnetic fields — Sun: photosphere — Sun: chromosphere 



1. INTRODUCTION 

The solar chromosphere is sandwiched between the sur- 
face, or photosphere, and the hot and tenuous outer 
corona. This highly structured region, on average 
2000 km thick, is constantly perturbed by short lived 
(3 - 10 minutes), jet-like extrusions that reach heights of 
2000 - 10000 km above the photosphere. These thin jets 
are formed in the vicinity of photospheric magnetic field 
concentrations. Until recently, their small size and short 
lifetinies have made detaile d analysis difficult ((BeckersI 
[ilHl iSuematsu et al.||199.^, which has led to a multi- 
tude of poorly constrained theories of their formation 
(fSterling 120001) . In addition, there has been consider- 
able confusion about the relationship between spicules 
at the quiet Sun limb, mottles observed on the quiet 
Sun disk, and dynamic f ibrils (DFs) found in the vicin- 
ity of active region plage (|Grossmann-Doerth Schmidtl 
11992*) , although the similarity in many of their properties 
strongly suggests some of these phenomena are related 
(jTsiroDoula et al.l 1199 4). We focus on observations of 
DFs (§2), compare them to advanced numerical simula- 
tions (§3), report on regional differences of DP properties 
(§4), and finish with a comparison to quiet Sun jets (§5). 
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2. OBSERVATIONS OF DYNAMIC FIBRILS 

The rece nt advent of the Sw edish 1-m Solar Tele- 
scope (SST, 'Scharmer et al.""2003l and advances in post- 
processing techniques ( van Noort et al. 2005) have al- 
lowed us to obtain an unprecedented, diffraction-limited 
(120 km) 78 minute long time series of the chromosphere 
as imaged in the core of Ha (656.3 nm.) at a cadence of 1 s 
(|van Noort Rouppe van der VQQrtll200 6^. These data, 
taken on 4 October 2005, resolve for the first time the 
spatial and temporal evolution of DFs, in particular the 
properties of the 257 DFs chosen for this study. While 
DFs have varying lifetimes, lengths and widths, a typical 
fibril rises rapidly to a maximum length in 1.5-3 minutes 
and recedes in a similar time along the same, relatively 
straight path, presumably parallel to the direction of the 
magnetic field (Fig. 1, top panels). Typical fibril life- 
times are between 120 and 650 s, with an average of 290 
s. DFs display some internal structure (e.g. at t = 139 s 
and t = 186 s in top panels of Fig. 1): many DFs do not 
rise and fall as a rigid body but rather show phase and 
amplitude variations in velocity at various positions away 
from the fibril axis. Despite the substructure, fibrils are 
thin, with widths ranging from the diffraction limit of 
120 km to 700 km. The maximum projected extent is 
usually relatively modest, ranging from 400 km to 5,000 
km with an average of 1,250 km. Both the lifetimes and 
lengths of DFs are in the lower r ange of values reported 
for quiet Sun spicules or mottles (|Beckers|ll968t) . 
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Fig. 1. — Temporal evolution of dynamic fibrils from Ho; linecen- 
ter observations at the Swedish 1-m Solar Telescope (top panels), 
and from numerical simulations (bottom panels). The observations 
show a dark elongated feature with an upper chromospheric tem- 
perature of less than 10,000 K rise and fall within 4 minutes. The 
bottom panels show the logarithm of the plasma temperature T, 
set to saturate at log T = 4.5, from numerical simulations covering 
the upper convection zone (z<0) up through the corona (white re- 
gion at the top). The horizontal scale is arcsec in the top panels, 
and Mm in the bottom panels. The vertical scale has its origin 
at the photosphere (optical depth rsoo = !)• Contours of plasma 
/3 are drawn in white where ^ = 0.1, 1 (thicker), 10. The plasma 
(3 is the ratio of the gas pressure to the pressure exerted by the 
magnetic field. 



It has been a subject of significant debate whether 
quiet Sun jets such as spicules or mottles follow paths 
that are parabolic, ballis tic (i.e., s ol ar gravi ty), or 
constant velocity (e.g ., iBeckersI fl968t tNishikawal 119881: 
ISuematsu et al .ill 9951) . The high cadence and high spa- 
tial resolution of our observations reveals that most ac- 
tive region DFs observed in Ha line core follow almost 
perfect parabolic paths. We have therefore fit the paths 
to all 257 DFs outlined in Ha linecenter with parabolas, 
and use these fits to determine velocities and deceler- 
ations. DF motions may be summarized by an initial 
impulsive acceleration, after which the top of DFs are 
subjected to a constant deceleration throughout their 
lifetime (Fig. 2, top panel). At the beginning of the 
ascending phase, the velocity of the fibril top is super- 
sonic, with an average (line-of-sight projected) value of 
18 km/s and a range from 10 to 35 km/s for the 257 DFs. 
The velocity of the fibril top then decreases linearly with 
time until it reaches a maximum downward velocity that 
is roughly equal in amplitude to the initial upward ve- 
locity. The average projected deceleration is 73 m/s^ 
with a significant spread, ranging from 20 to 160 m/s^. 
These projected values are significantly smaller than the 
(downward) solar gravitational acceleration of 274 m/s^. 
In principle, the deceleration suffered by DFs along their 
path could be much larger than the projected decelera- 
tion, if their path is close to the line of sight. We find 
that the projection angles necessary to obtain solar grav- 
ity (or its component along the field) are not compatible 
with the visual appearance of the large scale topology of 
the active region (Fig. 4, top panels). 

3. NUMERICAL SIMULATIONS OF DYNAMIC FIBRILS 

We find features with similar properties in recent ad- 
vanced two-dimensional numerical simulations (Fig. 1, 
bottom panels). These simulations model the evolu- 
tion of a radiative MHD plasma, and, for the first time. 



Fig. 2. — Space-time plots of: the height of several DFs from Ha 
linecenter observations (top panel), the logarithm of the plasma 
temperature in numerical simulations (middle panel), and simulta- 
neous plasma velocity from the same simulations (bottom panel). 
Most of the DFs in the top panel follow a near perfect parabolic 
path, as illustrated for one DF by the two separate fits (full and 
dashed lines) used to derive parabolic parameters. Parabolic paths 
with similar parameters are traced by fibril-like features in the sim- 
ulations (middle panel). The vertical scale, plasma (3 contours and 
temperature range are the same as in the bottom panels of Fig. 1. 
Upward plasma velocities (blue in the bottom panel) show the up- 
ward propagation of the shocks that drive the fibril-like features. 



span the entire solar atmosphere from the upper con- 
vection zone to the lower corona, including t he photo- 
spher e, chromosphere and transition region fHans teenl 
l2005l ). This model includes non-grey, non-local ther- 
modynamic equilibrium radiative t ransport in the pho- 
tosph ere and chromosphere ( N ordlundl Il982t iSkartlienI 
120001) . and optically thin radiative losses as well as mag- 
netic field-aligned heat conduction in the transition re- 
gion and corona. We used an average field strength of 
about 100 G, with the field clumping up to 1 kG in down- 
draft regions in the photosphere. In these models, acous- 
tic waves generated in or near the photosphere are found 
to propagate upward along the magnetic field lines, form 
shocks in the middle chromosphere, and lift the upper 
chromosphere several thousand kilometers along a front 
with a width of up to 1000 km. The elevated plasma pro- 
trudes into the hot corona (Fig. 1, bottom panels) and 
resembles the observed DFs (Fig. 1, top panels). The 
duration of these simulated jets is of order 3-5 minutes, 
with maximum lengths ranging from 500 km to several 
thousand km. 

Our numerical simulations explain why the DFs in our 
SST data have parabolic paths with decelerations that 
are only a fraction of solar gravity. In the simulations, 
fibril-like features are seen to follow parabolic paths with 
properties very similar to those of the observed DFs 
(Fig. 2, middle panel). These paths occur as a natu- 
ral result of shock wave driving. When a shock impacts 
the top of the chromosphere, the plasma is catapulted 
upward at a velocity that exceeds the local sound speed 
of ~10 km/s. This agrees well with the lower cutoff at 10 
km/s of maximum velocities in the observed DFs. Shock 
waves generally have velocity profiles in the form of 'N'- 
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or 'sawtooth' shapes (|Mihalas fc Mihalail 19841) . As a re- 
sult, a plasma parcel passing through a shock wave will 
first experience a sudden impulse in velocity, followed by 
a gradual, linear deceleration as the shock recedes — this 
is consistent with a parabolic path. 

The deceleration in these shock waves depends on the 
component of solar gravity that is parallel to the mag- 
netic field lines along which the shock waves propagate. 
However, it also critically depends on the period and 
amplitude of the shock waves. Even for vertical propa- 
gation, our simulations show that the deceleration is less 
than solar gravity. 'N'-shaped shock waves propagating 
at any angle to the vertical, will, for a given period, show 
a steeper decline with time (i.e., deceleration) in velocity 
for a greater shock strength. Similarly, for a given ampli- 
tude, shock waves with shorter periods show larger de- 
celeration. These dependencies explain why the observed 
DFs have decelerations less than solar gravity. In addi- 
tion, they offer a natural explanation for an intriguing 
linear relationship between the deceleration and maxi- 
mum velocity of the DFs in our SST data. We find that 
DFs with a larger deceleration show a larger maximum 
velocity (Fig. 3, top panel). This previously unknown 
correlation is well reproduced by the jets in our simu- 
lations (Fig. 3, lower panel), for which we find a linear 
relationship with a similar slope and similar range in val- 
ues as for the observed DFs. The observed relationship 
and the close similarities to simulations strongly suggest 
that DFs are driven by chromospheric shock waves with 
strengths between Mach 1 and Mach 4. 

4. REGIONAL DIFFERENCES IN FIBRIL PROPERTIES 

This scenario is further strengthened by the clear re- 
gional differences of dynamic fibril properties. The field 
of view of the SST observations contains a dense plage 
region (region 1, red circle in Fig. 4), where the mag- 
netic field is more vertically oriented, and a less dense 
plage region (region 2, blue circle in Fig. 4), where the 
field is more inclined away from the vertical. The dense 
plage region has a prevalence of DFs with larger decelera- 
tions (typically 100 m/s^, lower left panel of Fig. 4) than 
the region of less dense plage where DFs have smaller 
decelerations (typically 50 m/s^). In addition, DFs are 
shorter, typically 1000 km, and have shorter durations 
(3 minutes) in the dense plage region than in region 2, 
where DFs are on average 2000 km and typically last 
for 5 minutes. Since the regional differences imply high 
decelerations for short-lived and short fibrils, projection 
effects cannot explain these differences. If we assume 
that the low decelerations in region 1 are caused by a 
different viewing angle than in region 2, this projection 
effect would make the already longer fibrils in region 1 
even longer than those in region 2, which would make the 
regional differences more pronounced. Moreover, differ- 
ences in duration are independent of projection effects, as 
are various correlations between deceleration, maximum 
length and velocity. 

We believe instead that the regional differences occur 
because the solar atmosphere produces chromospheric 
shocks of varying strengths and varying periods. The 
strengths and periods depend critically on which dis- 
turbances propagate from the photosphere into the up- 
per atmosphere. The chromosphere acts to filter out 
upward-propagating disturbances with periods that are 
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Fig. 3. — Top panel shows for the 257 DFs observed at the 
SST the maximum velocity versus deceleration, both corrected for 
line-of-sight projection, assuming that the fibril is aligned with the 
local magnetic field as deduced from potential field calculations. 
These corrections are uncertain, but as the same correction is ap- 
plied to both quantities the correlation between them should not 
be affected by any errors this correction introduces. DFs in region 
1 and 2 are indicated by, respectively, red stars and blue diamonds 
(see Fig. 4 for definition of regions). The same scatterplot (lower 
panel), based on analysis of fibril-like features in the numerical 
simulations, reveals that the simulations reproduce the observed 
correlation between these parameters, as well as reproducing the 
range in deceleration and maximum velocity. 



longer than the local acoustic cutoff period. In gen- 
eral, the acoustic cutoff period depends on t he inclina- 
tion of the magnetic field line s to the vertical (|Suematsul 
119901: iDe Pontieu et al.ll2004D . Under conditions of ver- 
tical magnetic field such as in the dense plage region 
(region 1 in Fig. 4), this implies that the chromosphere 
is dominated by oscillations and waves with periods at 
the acoustic cutoff of 3 minutes. Since these waves shock 
and drive the DFs, this explains why DF lifetimes are 
about 3 minutes in region 1. Wavelet analysis of the 
Ha linecenter time series confirms that the upper chro- 
mosphere is dominated by wave trains with periods of 
order 3 minutes (Fig. 4, lower right panel). Observa- 
tions of the DFs in this region also show higher decel- 
erations and lower velocities. Fibrils in the dense plage 
region are driven by shocks with slightly lower amplitude, 
since the photospheric power spectrum peaks at periods 
of 5 minutes, and in addition the strong magnetic fields 
in dense plage regions generally reduce the amplitude 
of convective fiows and global oscillations at the photo- 
spheric level. These DFs also experience the component 
of gravity along the magnetic field, which in dense plage 
regions with more vertical field implies higher deceler- 
ations. The short lifetimes, high decelerations and low 
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Fig. 4.— SST images taken in Ho; wideband (0.8 nm FWHM, 
upper left) and Ho; linecenter (upper right) of part of NOAA AR 
10813, at heliocentric coordinates S7, E37. Top panels illustrate 
solar north and west, and the line-of-sight vector (0 = 39°). Blue 
and red contours in the upper left panel outline positive and nega- 
tive magnetic flux (± 110 Mx/cm^ in full-disk MDI magnetogram). 
Minor tickmarks in arcseconds. The DFs are predominantly ob- 
served in the unipolar plage region to the lower right of the sunspot. 
The longer, more static and heavily inclined flbrils, visible at, e.g., 
x=30" and y=45" are not the subject of this paper. The lower 
left panel shows an Ho; image superposed with dots, color-coded to 
illustrate flbril deceleration, at the location of each of the 257 DFs. 
The lower right panel shows the wave period containing the high- 
est number of wavepackets with signiflcant power in Ho; linecenter. 
Circles in the botttom panels show regions containing strong plage 
(region 1, red) and plage with inclined fleld (region 2, blue). 



velocities lead to shorter DFs. 

In region 2, where the field is inclined from the vertical, 
the acoustic cutoff period increases to 5 minutes, allowing 
waves with the full photospheric peak power at 5 minutes 
to propagate along the field into the chromosphere, and 
develop into shocks that are stronger than in region 1 and 
which drive DFs with lifetimes of 5 minutes. Wavelet 
analysis of region 2 confirms a preponderance of wave 
trains with periods of ~ 5 minutes (Fig. 4, lower right 
panel). The longer lifetimes and less vertical field lead 
to the observed lower decelerations and longer lengths. 

5. DISCUSSION AND CONCLUSION 

There are many striking similarities between quiet 
Sun mottles and the active region DFs studied here. 
Both phenomena appear as highly dynamic, dark fea- 
tures in the wings and core of Ha, and are associated 
with mag:netic fl ux co ncentrations. More importantly, 
iSuematsu et al.l (|1995D found evidence that quiet Sun 
mottles also follow parabolic paths with decelerations 
that are too small to be consistent with a purely bal- 
listic flight at solar gravity. While the interpretation of 
their observatio ns proved difficult wit hout detailed nu- 
merical models, S uematsu et al.l (|1995j) note that the ap- 
parent velocity profiles in mottles are fully compatible 



with impulsive acceleration followed by a constant decel- 
eration, with maximum upward velocities usually about 
equal in amplitude to the maximum downward veloc- 
ities. The velocities they report for mottles, of order 
10-30 km/s, are si milar to those we find i n our SST ob- 
servations of DFs. ISuematsu et al.l (|1995D also find that 
the largest Doppler velocities in mottles appear at the 
beginning of the ascending phase (blue-shifts) and at the 
end of receding phase (red-shifts), with downward red- 
shifted motion sometimes occurring close to their base 
during the ascending phase. These observations agree 
well with the properties of our simulated jets (Fi^. 2 , 
lower panel). In addition, iChristopoulou et al.l (|2001f) 
observe limb spicules with clear parabolic paths with de- 
celerations and maximum velocities similar to those for 
mottles and DFs. 

All of these strong similarities between previous mot- 
tle and spicule observations and our modeling and SST 
observations of DFs seem to imply that highly dy- 
namic chromospheric shock waves cause significant up- 
and downward excursions of the upper chromosphere 
in both active re gion and quiet Sun, as proposed by 
iDe Pontieu et al.l ((2004). Some unresolved issues remain, 
such as the longer lifetimes of quiet sun mottles and 
spicules (2-10 minutes), and the greater heights of 2-10 
Mm that spicules reach at the limb. Preliminary analysis 
of our simulations suggests that these differences could 
be related to large scale differences in magnetic topol- 
ogy. Further numerical simulations of various magnetic 
topologies will help resolve these issues. For example, 
it is possible that spicules reach slightly greater heights 
because they consist of two populations: jets that are 
driven by shocks (as described here), and jets caused by 
reconnection. The latter jets could form a subset that 
on average is taller than the shock driven jets, and per- 
haps be part of a continuous spectrum of reconnection 
jets th at includes surg:es , macrospicules and Ra upflow 
events (|Chae et al.lll998D . Whatever the role of recon- 
nection in quiet Sun, our findings indicate that, at least 
in active regions, most jets are caused by chromospheric 
shocks driven by convective flows and oscillations in the 
photosphere. 
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